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We examine the asymmetric simple exclusion process with open boundaries, a paradigm of driven 
diffusive systems, having a nonequilibrium steady state transition. We provide a full derivation 
and expanded discussion and digression on results previously reported briefly in M. Depken and 
R. Stinchcombe, Phys. Rev. Lett. 93, 040602, (2004). In particular we derive an exact form for 
the joint probability function for the bulk density and current, both for finite systems, and also in 
the thermodynamic limit. The resulting distribution is non-Gaussian, and while the fluctuations 
in the current are continuous at the continuous phase transitions, the density fluctuations are dis- 
continuous. The derivations are done by using the standard operator algebraic techniques, and by 
introducing a modified version of the original operator algebra. As a byproduct of these consid- 
erations we also arrive at a novel and very simple way of calculating the normalization constant 
appearing in the standard treatment with the operator algebra. Like the partition function in equi- 
librium systems, this normalization constant is shown to completely characterize the fluctuations, 
albeit in a very different manner. 
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I. INTRODUCTION 

Since it is typically not feasible to describe macroscopic 
complex systems in terms of the exact dynamical evolu- 
tion of their microscopic degrees of freedom, our knowl- 
edge of the collective properties of such systems has come 
mainly through using the methods of statistical mechan- 
ics. The systems are there described by a reduced set of 
variables (not necessarily in a direct correspondence to 
any physical observables) evolving according to stochas- 
tic dynamics, modeling the effect of the suppressed de- 
grees of freedom. For equilibrium systems this approach 
has been extremely fruitful. On the other hand, recent 
insights into the properties of non-equilibrium steady- 
states (NESS) have come largely from specific micro- 
scopic studies (see below) and our general understand- 
ing of these steady states is still far inferior to that of 
equilibrium steady-states (ESS). 

On a phenomenological level, the distinction between 
ESS and NESS lies in that only the latter allows for net 
probability currents through the state-space (in equilib- 
rium detailed balance excludes this). The currents are 
set up and maintained by contacts with multiple reser- 
voirs at different potentials (e.g. different temperature 
or chemical potential), and/or by the presence of non- 
conservative bulk forces. We will refer to the former case 
as a boundary driven system, and to the latter as a bulk 
driven system. On a formal level the difference between 
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a ESS and a NESS lies in that for ESS the probabil- 
ity weights of different configurations are governed by 
detailed balance, while there is no such generic rule for 
NESS. Consequently the probabilistic weights of the con- 
figurations in a NESS are not a. priori, known (contrast 
the Gibbs weights of ESS), and the lack of detailed bal- 
ance allows for net-currents through the configuration 
space. Apart from this, there is also the fact that ESS 
with short ranged interactions generically (i.e. in the 
absence of a spontaneously broken continuous symme- 
try, and not at a continuous phase transition) have finite 
correlation lengths. Therefor one can adopt an ensemble 
approach to fluctuations for many-particle ESS [l| . For a 
NESS on the other hand, there can exist long-range cor- 
relations even when the system is not externally tuned to 
be at a phase transition jyj . This observation has lead to 
the ideas of self-organized criticality, whereby the system 
generates infinite correlation lengths by self-adjusting its' 
state to a critical state 0- By now there are also exam- 
ples of non-critical NESS, which still show infinite corre- 
lation lengths [J [E 01 or (as in the system studied below) 
where boundary driving is felt throughout the system . 

Faced with these complications one does wisely in uti- 
lizing methods that have proved their worth elsewhere. 
One highly effective tool within the study of ESS is the 
use of specific simplified and abstracted models. Among 
these are the classical- and quantum-magnetic models, 
the Ising model, and the Heisenberg models to name 
only two. These models were initially only hoped to give 
a qualitative description of the physical systems. They 
retain only some of the essential ingredients, while leav- 
ing out enough of the details for progress to be possible. 
With the advent of renormalization group arguments , 
and the implied universality, these models are now known 
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to also give quantitative predictions, and the study of 
simple models is now a fundamental part of the under- 
standing of ESS. 

Led by this, and the fact that universality seems to 
be a generic feature also of the long wavelength and low 
frequency degrees of freedom of NESS much of the ef- 
fort within the study on NESS has focused on particular 
simple models. Here we study a well know simple 1 + ID 
model of a driven diffusive system, the asymmetric sim- 
ple exclusion process (ASEP). It belongs to the class of 
stochastic interacting particles models [3, and the par- 
ticular version we choose to study is driven both from 
the boundary and the bulk. This model is non-trivial, 
displaying macroscopic collective behavior in the form 
of both continuous and discontinuous boundary-induced 
phase transitions. It is yet simple enough to be inte- 
grate pjj m mm. It is of additional interest since 

it maps onto certain growth models 0, it models traffic 
flow [l4[ , and it is believed to describe the la rge scale dy- 
namics of the noisy burgers equation \13L \WL fl7| and the 
KPZ equation [lEG3. There has been much progress 
in the analytical treatment of the ASEP, giving rise to a 
host of exact results describing its' steady state proper- 
ties [ISIlllllEilllllJIIJil. The probability mea- 
sure of this model (and of certain other one-dimensional 
models (see e.g. section 6 of [lTj)) turns out to be a 
generalization of the trivial product measure, where in- 
stead of a measure built up of commuting c-numbers at 
each spatial position, the measure is constructed by as- 
signing non-commuting operators to each spatial position 
(for further details see below) . Even though the fact that 
the system maintains a current is ultimately what sets it 
apart from an equilibrium system, the results so far con- 
cerning the currents are mainly for systems with periodic 
boundaries or infinite geometries with special initial con- 
ditions mm mm. 



In this paper we expand the discussion and digress on 
what was briefly reported in (25|. We will derive the 
exact joint probability distribution for the system av- 
eraged current and density, as well as their asymptotic 
form in the thermodynamic limit. This is a step closer 
to the achievements of statistical mechanics for ESS, and 
comparable to summing the partition function; but it is 
system-specific. 

The paper is organized as follows; Section ITT1 gives the 
definition of the model, recapitulates a few known results, 
as well as giving a brief introduction to the operator al- 
gebraic approach. This is followed by Section ITU where 
we introduce a novel way of calculating the normalization 
factor (first introduced in Section IIIBfl , and further in- 
troduce a "relaxed" operator algebra, with which help we 
are able to calculate the exact form of the density-current 
probability function. This is complemented by the study 
of the thermodynamic limit. In Section llVI we translate 
the results into the language of the KPZ-equation. The 
details of the calculations of the preceding sections are 
recorded in the appendices. 



II. MODEL DEFINITION AND SOME KNOWN 
RESULTS 

We consider the totally asymmetric exclusion process 
(ASEP) on a finite chain of size L with open boundaries. 
The site label I runs from left to right (see Figure ^1 . 
Each site on the lattice can be occupied by no more than 
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FIG. 1: Illustration of the dynamical rules of the ASEP 

one particle, n/ € {0, 1}. Given that the right neighbor- 
ing site of an occupied site is empty, the occupying par- 
ticle will jump to the empty site with rate 1. If the first 
site on the lattice is unoccupied, particles are injected at 
this boundary with rate a. Further given that we have 
a particle at the last site of the lattice, it is ejected with 
the probability rate 0. No further transitions are allowed. 
We will here limit our considerations to the case where 
we can view the boundary rates as deriving from parti- 
cle reservoirs. We therefore take < a = pi e ft < 1 and 
< = 1 - flight < 1, where p\ c{t and p right are the 
particle densities of the reservoirs. 



A. The phase diagram 

This model has been exactly solved ^1 m the sense 
that the steady-state probability of any given microscopic 
configuration can (in principle) be calculated by apply- 
ing a given set of algebraic rules (described below) . This 
though is far from what we are used to refer to as a 
solution in equilibrium statistical mechanics. There we 
are in a better position already from the start in that 
the configurational weights are explicitly given, while for 
this model they are given only in terms of algebraic rules 
(or matrix multiplication in case we have a finite ma- 
trix representation of the operator algebra (see below)). 
In equilibrium it is normally the successful summation 
of the partition function that is considered a solution 
of the problem. Even so, the algebraic rules yield a 
wealth of information about the system. Especially they 
can be used to deduce the precise form of the phase- 
diagram If consists of three parts A, B, and C, 
as given in Figure [21 In part A, the low-current, high- 
density phase, the average bulk profile, pi = (n/), is flat 
in the bulk with pi = p r ight = 1 — 0, and the current 
is ji = 0(1 — 0). The average bulk profile is completely 
dictated by the right hand side reservoir. The bulk pro- 
file connects to the value of the left reservoir through an 
exponential decay, with some characteristic length. The 
situation is reversed in B, the low-current, low-density 
phase, in that the bulk density profile is dictated by the 
left reservoir, with p; = pi c f t = a, and ji = a(l — a). 
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point is to represent any microscopic configurations in 
terms of a string of non-commuting operators D and E, 
corresponding to a particle and a hole respectively. It can 
then be shown that the steady-state probability function 
can be written in terms of this operator string and two 
auxiliary vectors, (a\ and \0), according to 



P ss ({ ni }) = {Zf r l {a\X{ ni )X(n 2 ) 



■X(n L )\f3). (1) 



Here the operator X(ni) equals D if there is a particle at 
site I (ni = 1), and E if site I is unoccupied (ni = 0). The 
vectors (a\ and \0) describe the properties of an uncor- 
rected particle reservoir. The state independent factor 
Z^f = (a\(D + E) L \(3) ensures the proper normalization. 
For to hold true, the operators and vectors must fur- 
ther satisfy the algebraic rules 



FIG. 2: Phase diagram of the one dimensional exclusion pro- 
cess. The dashed and dash-dotted lines indicate respectively 
first order and continuous transition lines. 



The decay of the density profile to the right boundary 
value /9 r ight = 1 — /? is also here exponential. Separating 
the two phases is a phase transition at which the typi- 
cal bulk profile develops a kink, taking the density from 
the value of the left reservoir to the value of the right 
one. The kink is of finite extension, and equally likely to 
be situated anywhere throughout the bulk. When pass- 
ing over this transition line the density is discontinuous, 
and the border between A and B is thus a first-order 
non-equilibrium phase-transition. Such a boundary in- 
duce phase transition has no counterpart in equilibrium 
systems where boundaries are assumed irrelevant. The 
remaining region, region C, is called the maximal-current 
phase. Here the bulk properties are set by the bulk drive, 
i.e. in this phase both the injection and ejection rate, are 
high enough for the system only to be limited by the 
transition rates in the bulk. The bulk profile is given by 
pi = 1/2 and ji — 1/4. In this phase the decays from the 
boundaries are algebraic. 

A nice heuristic argument for the form of the boundary 
decays has recently been presented in [2(|. The transi- 
tions from A or B to C induces a jump in neither the av- 
erage density nor current, and thus these transitions are 
continuous. As we will argue later, they are in fact second 
order. The point where all regions meet is called the criti- 
cal point. These results can all be derived by utilizing the 
matrix algebra outlined below. Further results contain- 
ing information about the time evolution of the system 
can be derived by Bethe-ansatz methods [T(| |2(J, 123, £3 • 



B. Operator algebra 

We here outline the operator algebra as given in [ll| 
since we will be using and generalising this. The starting 



DE = E + D d = C, {a\E = -{a\, D|/3) = ±|/3>, (2) 

a p 

and we may take the normalizations of the vectors \a) 
and (j3\ to be such that (a\j3) = 1. The algebraic rules (j2J) 
are now all that is needed to calculate P^{{ni}), result- 
ing in a polynomial of degree L in 1/a and 1//3. An al- 
ternative approach is to look for matrix representations 
of the algebra ©, and then use these to calculate the mi- 
croscopic probabilities. This justifies calling these states 
matrix-product states. Even though we now have a way 
of getting the steady-state weight of any specific configu- 
ration, actually calculating this number becomes increas- 
ingly hard as one considers larger and larger systems. 
Thus one wishes to extract general information directly 
from the algebraic rules, without explicitly calculating 
the microscopic weights. This is done for both average 
density and current in ^lj, and we here just recall that 
it is in general very easy to write down the desired quan- 
tities in terms of the above defined operators 

Pl d ^{n l ) = (a\C l - l DC L - l \(3)/zf 

ji = ( n i0- — n i+i)) 

= {a\C l - 1 DEC L - l - 1 \f3)/Z2 = Zf_jzf 
C(l,m) d = (nin m ) 

= (a\C l - 1 DC m - l - 1 DC L - m \p) /Zf . 

The first two can be calculated asymptotically for large 
L and fixed x = n/L [Tl| . giving the phase diagram dis- 
cussed above. Many of the results for these systems have 
been derived through finding a matrix representation of 
the operator algebra. If there exists a finite sized matrix 
representation, then the above form of the correlation 
functions show that the inverse correlation length is sim- 
ply proportional to the highest eigenvalue of the matrix 
representation of C. Thus, only when no finite repre- 
sentation can be found is it possible to have an infinite 
correlation length, and algebraic decay of correlations. 
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III. MACROSCOPIC DESCRIPTION 

As is the case for ESS, the transition from a detailed 
microscopic knowledge about the weights of each con- 
figuration to information concerning macroscopic prop- 
erties of the system is in general highly non-trivial (c.f. 
summing the partition function). Here though we are 
faced with one additional problem in that we do not 
have the microscopic weights explicitly, but only a set 
of algebraic rules J5J), or as a string of possibly infinite 
matrices. Thus, as mentioned above, we need to extract 
macroscopic quantities directly through using the alge- 
braic rules. One macroscopic quantity that has been 
likened to a partition function of this non-equilibrium 
system is the normalizing constant Z^f |29j . Though it 
lacks the power of yielding moments through differenti- 
ation, it generates a form of Lee- Yang theory of these 
non-equilibrium phase transitions. Through the consid- 
erations in this section we will further see that the nor- 
malizing constant still plays an instrumental role in de- 
termining the complete statistics of the bulk densities 
and currents in this non-equilibrium system. 



A. The generating function for 

As a warm up to what will follow, and since the nor- 
malization constant plays a central role in later develop- 
ments, we here present a novel and very simple way of 
calculating Z^f . This is done through calculating the 
generating functional of the normalization constant, an 
approach recently used in |3(], On] . We define the gener- 
ating functional as 
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where for an arbitrary operator X we have defined 

x 
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(4) 



Here X is some finite integer, and it immediately follows 
that 



(1-fiX) 



1 



= l + 0(p x+1 ) 



x 



1-(J,X 

From the operator algebra J2J we have 

(i-nD)(i-»E) = i-/z(i-m)c, 

and inverting this in the sense of JHJ we get 



(5) 
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1 - fiE_ 
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1 - fiD 



-0(/x 



A + lN 



A 



In the above we note that all .E's are to the left of all D's. 
Since (a\ and \[3) by definition are respective eigenvectors 



of these operators, we have 
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1 — [ij a 
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[l-fx/(3\ 



+ 0(M A+1 ) 
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In the limit X — > oo, the above expression has a conver- 
gence radius of min(a, 0) > 0. Thus as long as we are 
within this we can take this limit and write 



G^( M (1 - aO) 







a — /i (3 — n 



There now exists an open region around the origin where 
this can be rewritten as 



2a 



2/3 



2a - 1 + VI -4/i 2(3 - 1 + VI -4a*' 



The expression can be analytically continued to all //, 
and is inverted through 



Zf = 



d/i G al3 (p) 
2m n L+1 



(6) 



where encircles only the pole at the origin, and does 
so once in the positive direction. This can be used to 

a/3 



derive the same finite size form of Z 
l i 



7 aP 



i 



k Rl-k ' 



1=1 



k=0 



a k ft 



Ala = 



L as given m [Tl| . 
l{2L-l-l)\ 



L\(L-l)\ 



(7) 

One can also perform a large system asymptotic analysis 
on © using steepest-descent methods, which again yields 
the same results as in [TT| . 



1/2 <a< 13: zf 
a < 1/2, 0: Zf 
a = (3 < 1/2: Z 



A L+1 a(3 a + (3 -I 
L 3 / 2 V^ (2a- l) 2 (2/3- 1) : 



-,(8) 



/3(1 - 2a) 



1 



a/3 



{(3-a){l-a) a L {\-a) L 
(1 - 2a) 2 L 
(1-a) 2 a L (l-a) L ' 



(9) 



The asymptotic forms for a > [3 can be obtained through 
realizing that the system exhibits a particle hole symme- 
try. That is, instead of focusing on the particles we might 
just as well consider the holes as evolving with exactly 
the same dynamics, but with the injection and ejection 
rates exchanged. Thus we can directly get the result for 
a > f3 by letting a — > (3 and (3 — » a in the above. 



B. The bulk current-density probability function 

Though the normalization constant considered above 
has some of the feature of the equilibrium partition func- 
tion, it does not in it's present applications tell us much 
about the moments of the two natural observables of the 
system, the density and current. Thus we here concen- 
trate on the derivation of the exact joint probability func- 
tion for the average density and current throughout the 
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bulk. Through this we will see how the normalization 
constant also here tells us about fluctuations, albeit in a 
manner very different from that of equilibrium statistical 
mechanics. This is done for any system size, and later 
the thermodynamic limit is also considered. 

First we define the total activity within the system as 
the number of bulk bonds that can facilitate a transi- 
tion of a particle in the immediate future, i.e. the total 
effective bulk transition rate. The bulk current is then 
defined as the activity divided by the system size. For 
any given state the activity equals the number of pairs 
of neighboring sites that have a particle to the left and a 
hole to the right. To get a handle on the activity, J, of 
a microscopic configuration of N particles we choose to 
represent such a configuration by a sequence of J objects 
of the form D Pj E hj , Pj,hj > 1, possibly padded with 
E's to the left and D's to the right. Each of these ob- 
jects contains what corresponds to an active bond, and 



using these objects we can write any microscopic steady 
state measure as 

Pw({Pj,hj}) 

= (Z2 P y 1 (a\E ho {D pi E hl )- ■ ■ (D PJ E hj )D P0 \(3) 1 

by appropriately choosing the numbers {pj,hj} and J. 
It further follows that the above expression is unique if 
h-OiPo > Oj and the rest satisfy hj,Pj > 1. We can now in 
principle calculate the joint probability distribution for 
N and J by summing the above over all /ij's and p^s con- 
sistent with a specific number of particles (X) j=o Pi = -^0 

and a given system size (N + Y^,j=o ^i = L)- Choosing 
to enforce these constraints with contour integral repre- 
sentations of the Kronecker delta, the expression for the 
joint particle-activity probability function can be written 
as 



L N,L-N 

Pf(N,J)= £ J2 p ^{Pj^j})S Epj , N S Ehj , L . 

po : ho—0 pi,hi — l,l>l 



N 



a/3 



dzdi 



G Z ,C, 



(2m) 2 z N+l-JzL-N-J+l ( z _ ^ _ a ) 



1 - zD 



DE 



N-l 



1 - zE 



|/3). (10) 



L-N-l , 



Here C z (C 2 ) is a directed contour that encircle the pole 
at the origin of the complex z (z) plane once in the pos- 
itive direction, with \z\ < /3 (\z\ < a). The first step 
toward explicitly calculating (|10f) is through considering 
the properties of involved operator product. Surprisingly 
one can show (see Appendix^) that a slight modification 
of the above operators 



D' 
E' 



d ^[l-(z + z) 



del' 



[!-(* + *)] 



1-zD 
1 



1-zE 



D. 



N-l 



L-N-l 



satisfy the "relaxed" operator algebra 

D'E' = D' + E' + 0(z N ,z L - N 



(11) 



(12) 



The new relaxed eigenvectors and eigenvalues are simply 
given by 



D'\ 



1 



&(z N ), 



(a\E' = ± i (a\+0(z L - N ), 



with the eigenvalues defined as 



, dcf 



l-(z + z)' 







dcf 



l-(z + z)' 



(13) 



The fact that these eigenvalues are complex is of no con- 
cern since we consider only finite polynomials in the in- 
verse eigenvalues. Any result is thus uniquely extendable 



into the complex plane through analytic continuation. 
We can rewrite <|10[1 in terms of the primed operators, 
and start using the relaxed operator algebra to trans- 
form the expression. The result of any such manipulation 
would, according to the above, be the same up to terms 
of order z N and z L ~ N , as if the operator algebra would 
have been exact. Terms of this order have no effect under 
the contour integral in I|1U|) since the poles at the origins 
are of order equal to or lower than TV and L — N respec- 
tively. (The case for J = is trivial.) Thus, using the 
new algebra to perform any manipulation within (|l(Jfl is 
equivalent to using the exact algebra. Therefore we can 
write 



Pf(N,J) 



a(3 



dzdz 



1 



Cz,C s 



(2m) 2 z 1 



7 a'P' 



(z - (3){z - a)[l - (z + z)} 2J ' 



(14) 



This expression is the main result of this work, and since 
all quantities in it are known exactly, it yields both the 
exact finite system size form of P£ (N, J), as well as the 
asymptotic form in the large system size limit. Using 
the above expression it is further easy to derive a similar 
form for the joint generating functional of the density 
and current. This is outlined in Appendix iBl while we 
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go on here and present exact and asymptotic results for 
the probability function. 

C. Exact results for finite systems 

The integral in l|14fl is easily calculated with the help 
Cauchy's integral theorem. All we need to do is to cal- 
culate the coefficient of the term proportional to (zz) -1 
in the Laurent-series expansion of the integrand in (|14|l . 



For the special case J = we have Zj & = 1, and thus 
Pf(N,0) = ^(l//3f(l/a) L - N . 



This is obviously correct since the inactive state must 
have L — N empty sites followed by N filled sites. In 
Appendix |0 we consider the case J > 1. The result is 



„ J J L-N-JN-J 

L j=l k=0 c=0 d=0 



with the combinatorial factors 

'k + c\ 1 



G k , c ( a ) = 



c J a c+k+1 
K-l + a + e\fa + e 
a + e J \ e 



Through the above we now have the exact form of the 
sought-after joint probability function for any system 
size. The form is illustrated in Figure |3J 



D. Thermodynamic limit 

We here return to <|14|1 . Using the asymptotic form of 
the normalizing constant given in (jSJ) and ©, we perform 
a steepest-descent calculation to get the asymptotic re- 
sults in the large system limit. We consider the different 
phases individually. Due to the particle-hole symmetry 
P^(N, J) — P^ a (L — N,J), it is only necessary to ex- 
plicitly consider the case a < (3. 

First turning to the maximal-current phase we use i|14|) 
together with (JSJ) , and , and drop all pre- factors that 
are independent of N and J (this will be done through- 
out), to write 



Pf(N,J) 



dzdz 



1 



J 3 / 2 / {2m) 2 zN-J+l^L-N-J+l 

1 1 

X [1 - (z + z)] 2 ' 7 - 1 [2a - 1 + (z - z)] 2 [2/3 - 1 - (z ~ z)} 2 ' 

The asymptotic behavior of these integrals is in princi- 
ple straight forward to calculate. In practice though, it 
turns out to be quite cumbersome since one has to de- 
termine which of the saddle points and lower order poles 
give the dominant contributions. We can shortcut this 
through only considering the asymptotic form in some fi- 
nite region around the peak of the distribution. From the 



general discussions of the phase diagram in Section Til Al 
we know that the average density and current is a and (3 
independent. Thus, the lower order poles cannot dictate 
the asymptotic behavior around the peak value of the 
probability distribution, and instead this must be set by 
the saddle points 

z* = p - j, z* = l-p-j, p = N/L, j = J/L. 

A saddle-point approximation thus results in 



1 



(16) 



where we for simplicity have dropped all the sub- 
dominant pre-factors. Even though the extent of the re- 
gion of validity of (| 1 fcil) is unknown, it should be pointed 
out that the size of this region is a finite fraction of the 
complete range of p and j (as long as the system is away 
from any phase boundaries). In the first row of Figure 01 
we show the resulting dominating asymptotic plots. 

We now turn to the low current, low-density phase. 
Using |[nj, ©, and (0 we have 



Pf(P,3) 



dzdz 



, )0 , (2m) 2 z n-J+i e l-N-j+i 
2a - 1 + z- z 1 



/3 — a — (z — z) (a — z) J+1 (1 — z — a) 



J+i • 



The same arguments as applied in the high current phase 
gives the asymptotic probability distribution around the 
peak. Again it is the saddle points 



P ~ 3 \ -* 1 - P ~ 3 
z = (1 — a), z — — a, 



1-P 
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FIG. 3: (Color online) Each row contains a surface and a 
contour plot of the exact joint probability distribution for 
the values of a and f3 indicated, and with p — N/L and 
j = J/L. The first three rows illustrate the behavior of the 
probability distribution as the system goes along the line of 
a — [3 through the critical point at a = j3 = 0.5, while the 
last three graphs illustrate the behavior as the system goes 
through the first-order transition at a = j3 = 0.25. Overlaid 
in the contour plots (dashed line) is the curve j = p(l — p) 
which defines the set of possible asymptotic average values of 
p and j throughout the system's different phases (not at the 
first order transition line) . The system size is L = 40. 




FIG. 4: (Color online) The two rows display a surface and a 
contour plot of the leading behavior of the asymptotic joint 
probability distribution. The calculations were performed at 
the injection and ejection rates indicated and at a system of 
size L — 40 (to make the result comparable to Figure [3J. 



that dominate. The resulting dominant form is 

( p^l-p) 1 -" 1 V 

\a^P{l - a)P fi{p - j)P~i(l -p- jy-P-i ) 

(17) 

The above result is directly transferable to the high- 
density phase through the use of the particle hole symme- 
try mentioned above. A realization of the asymptotically 
dominating part in the low-density phase is shown in the 
second row of Figure 0] 

It is clear from the asymptotic forms that the prob- 
ability distribution is non-Gaussian in all phases. This 
is consistent with the view that long-range correlations 
are a generic feature of non-equilibrium systems with 
locally conserved dynamics 32] . Sufficiently close to a 
phase transition, any finite system will reach a point at 
which the region of validity of the above asymptotic forms 
shrink to the size of the typical fluctuations. When this 
happens the system crosses over to a situation where the 
fluctuations are governed by the tails excluded in the 
above development. 

It is interesting to note that if we consider the proba- 
bility density of the current at a fixed density, we recover 
the same functional form as derived in |3i| for the prob- 
ability density of the current in a closed periodic system. 

E. Fluctuations 

It is interesting to note that as the continuum tran- 
sition is passed, the asymptotic forms l|l(jD and (|17|l in- 
dicates that there will be a jump in the density fluctua- 
tions of the system. More precisely, in the low-current, 
low-density phase we have 




a(l - a) f 1 1 - 2a \ 

~ L \ 1 - 2a 1 - 3a(l - a) ) " 

Using the particle-hole symmetry we have 

for the low-current, high-density, phase. In the maximal- 
current phase we have 

r c If 1/8 \ 
mc ~ L \ 1/16 J " 

Thus as we go from either of the low-current phases (re- 
gions A or B in the phase diagram in Figure |2J) to the 
maximal current phase (region C in Figure |2J) , there will 
be a jump in the correlator 
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Since we here have a discontinuity in the density fluctua- 
tions as we pass over the continuous phase transition, we 
see that this transition is of second order (c.f. equilibrium 
statistical mechanics where a discontinuity in the corre- 
lator corresponds to a discontinuity in the second-order 
derivative of the free-energy). The sign of this jump also 
illustrates that the strength of the fluctuations decrease 
as we enter the maximal-current phase. 



IV. IMPLICATIONS ON THE KPZ EQUATION 

In this section we briefly point to the known connec- 
tions between the ASEP and the KPZ equation |sl Il7|. 
and translate our findings to the language of the KPZ 
equation. 

The ASEP can be mapped onto a lattice growth 
model 0, ^J, which in turn is believed to share its' 
long-wavelength characteristics with the KPZ equation. 
Within the framework of this standard mapping it fol- 
lows that the growth velocity and the average slope of 
the interface are respectively given by 

v = j, dh = 1 — 2p. 

Letting 9/i|i e ft /right denote the enforced boundary slopes, 
we can get the joint slope-velocity distribution for the 
interface model by substituting 

p = (i-aft)/2, j = v_ 

a = (1 - dh\ieit)/2, 13 = (l+M| right )/2 

into the probability distribution for the ASEP. It is fur- 
ther clear that when the average slope at the left bound- 
ary, dhieft = 1 — 2a, matches the slope at the right 
boundary, dh r i g ht — 2/3—1, i.e. when a + /3 = 1, we 
are at the trivial line of the ASEP where the matrix- 
product measure reduces to a product measure. Situa- 
tions where more complicated boundary conditions are 
relevant have recently been examined experimental and 
theoretical in EES. 



V. CONCLUSION 

In this paper we have examined the joint probability 
distribution of the system-averaged density and current. 
We have derived an exact expression for the joint prob- 
ability function for finite systems, and also considered 
the thermodynamic limit. This was done by introducing 
a relaxed operator algebra, and it would be very inter- 
esting to examine if the same "trick" could somehow be 
applied to the PASEP. This is especially important since 
this model interpolates between a ESS and NESS. The 
development further shows that even if the normalization 
constant does not act as a partition function in the nor- 
mal sense of giving moments through differentiation by 
a conjugate field, it nevertheless completely governs the 
fluctuations through 114(1. Sufficiently close to a phase 



transition, any finite system will reach a point where the 
above derived asymptotic forms are not valid. Thus it 
would be interesting to derive the full asymptotic form 
of the probability distribution, including the tails. We 
have also shown that the continuous transitions are sec- 
ond order in the sense of equilibrium statistical mechan- 
ics, and that it is the density fluctuations that display a 
discontinuity. Lastly we wrote down the translation of 
the probability density for the asymptotic ASEP to the 
probability density of the KPZ equation in terms of the 
average slope of the interface, and the average interface 
velocity. 
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APPENDIX A: THE RELAXED OPERATOR 
ALGEBRA 

The major problem with calculating the probability 
function as given in l(TU|) is the appearance of the product 
of the operators 



1 



1-zD 



D and 



A" 



1 



1-zE 



x 



where all the D's are to the left of all the E's. By con- 
sidering this product we will uncover operator relations 
very similar to those of the original operator algebra J3J) • 
Using the original algebra we can write, 



1 - zD 



DE 



1-zE 



zD 



x 
D 



1 



1 zD 



1-zE 
E 



x 



1 - zE 



(Al) 



x 



From the inversion relation Q we have 



1 



1 — /iX 



= 1 + 



x 



1 — fiX 



+ 0(^ x+1 ) 



x 



and using this in the right hand side of IjAlJ) it is easily 
seen that (|A1|) is equivalent to 



(1 



A) 



i 



1 - zD 
1 



1- zD 



D 



D 



1 



1-zE 
1 



x 



X 



1-zE 



x 



&{z 



X + l -X+l 



z A+1 ). (A2) 
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It is now clear that the operators (|llfl will satisfy the 
relation In terms of the operators the central 

operator product in I jlOjl reads 



We proceed by considering the expansion factor by factor 
in their respective Laurent-series around the origins. The 
first two factors are given by the expansion 



1 - zD 



DE 



N-l 



1 - zE 



L-N-l 



D'E' 



APPENDIX B: GENERATING FUNCTIONAL 



with 



G kiC (a)z c 



By using the definition of the generating functional 
for Z°lf (j2J, and the integral representation of the joint 
probability function Ijl4jl . we can write the generating 
functional for N and J as 



Gk, c (ot) 



k + c\ 1 



a(3 f dzdz I 



Zj P Jc z ,Cz ( 2m ) 2 ~ a )( z ~ I 3 ) z L+1 z-jz 



x G 



a'0 1 



{l-[z + z])< 



The contours should be chosen such that the only en- 
closed poles are the ones at the origins. This object has 
a rather complicated analytical structure and for sim- 
plicity we have restricted our efforts to the probability 
function. 



while the last factor expands as 



with 



(l-[z + z}) 



2J-1 



— E H2J-La,eZ e Z a , 
a,e— 



K — 1 + a + e\ fa + e 
a + e 



APPENDIX C: EXACT PROBABILITY 
FUNCTION 



Considering the product of all three factors we arrive at 



In order to get the exact probability distribution for iV 
and J we need to expand the integrand of (|14f) around 
the origins of both the z- and z-plane. This is done in 
the present section. First we write down the complete 
expression for the probability function, given the exact 
form of the normalizing constant as shown in Q 



a'0' 



Z 



Thus we have 



Pf{N,J) 



E^E 



k=0 



(!-[* + !])' 
{a-z) k {(i-z) l - k ' 



E^E 



J L 1=1 

X 



fc=0 



dzdz 



1 



a/3 

7^ Z_ ■'*-■< Jk (2m) 2 z n+1-J E l-n+i-j 



(a - z) fc +!(/3 - «)'-*+!(! - [z + z]) 2J - 1 ' 



(a - z) fc +!(/3 - (1 - [z + z]) 2J - 1 

CO 

= ^ G k,c( a ) G l-k,d(0) H 2J-l,a,, 
c,a,d,e—0 
oo oo 

= E E G k,c( a ) G l-k,d(P) H 2J-l,a- 
c,d—0 a—c^e—d 

oo / a.e 

= ^2 I ^2 ^k,c(^)Gi-k i d(P)H2J-l >a -c,e- 
a.e=0 \c,d=Q 



-c+a z d+e 



c,e—d% % 



d I z a z e 



The only terms in this series that will contribute to the 
probability density are the ones with a = N — J and 
e = L — N — J, and thus we have the exact form of the 
probability function given by JTHJl. 
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